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PREDICTION OF COMPONENTS IN RANDOM SUMS 

MUNEYA MATSUI 

Abstract. We consider predictions of the random number and the magnitude of each iid component 
in a random sum based on its distributional structure, where only a total value of the sum is available 
and where iid random components are non-negative. The problem is motivated by prediction prob¬ 
lems in a Poisson shot noise process. In the context, although conditional moments are best possible 
predictors under the mean square error, only a few special cases have been investigated because of 
numerical difficulties. We replace the prediction problem of the process with that of a random sum, 
which is more general, and establish effective numerical procedures. The methods are based on con¬ 
ditional technique together with the Panjer recursion and the Fourier transform. In view of numerical 
experiments, procedures work reasonably. An application in the compound mixed Poisson process 
is also suggested. 


1. Preliminaries 

Motivated by prediction problems in a Poisson shot noise process, we consider two types of 
problems for random sums of iid random variables (r.v. or r.v.’s for short). Let A be a non¬ 
negative integer-valued r.v. and denote an iid sequence of non-negative r.v’s by 2,... so that 
SN = denotes the total sum. The distributions of both N and Xi are assumed to be known. 

Our problem is how we could obtain the information of the number N or each component A, 
when we only observe S n. Although there are several methods for these quantities such as linear 
predictions cS with c some constant, our methods are those by conditional moments, which are 
minimizors of the mean square error. More precisely our focus is on the following two types of 
conditional moments: 

(1.1) and | S/v], forkeN, 

where (Xi) may take both real and integer values and N denotes the set of natural numbers as usual. 

This type of random sum S ^ has been studied for a long time and has applications in a variety 
of fields. One could find many examples in the book of Feller |I5] XII] such as genetics, required 
service time, cosmic ray showers, and automobile accidents to name just a few. A large num¬ 
ber of relevant researches have been conducted, including e.g. calculations for probability ofS^ 
(Sundt and Vemic |[28]l ) or various limit theorems (see e.g. Gut |[8l and consult a nice summary in 
Embrechts et al. dH 2.5]). In recent years tail asymptotics have intensively studied, since accurate 
calculations of tail probabilities of S are computationally quite expensive, while they are required 
in applications. See lessen and Mikosch [|9l| for a survey with regularly varying tails and Goldie 
and Kliipperberg @ for that with subexponential tails. 

In this paper we do not go further into asymptotics but investigate precise calculations of quanti¬ 
ties (11.11) . which have not been studied yet except for some special cases (see Subsection ll.il) . We 
rely on two numerical methods, i.e. the Panjer recursion and the Fourier method, which are useful 
tools for computing P(Sn = n) and which are competitive ([|3||). The Panjer recursion scheme 
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originated in Panjer Il2^ is known to be stable when N belongs to the Panjer elass in most oases 
(lUSII). Meanwhile, the Fourier method oould be applieable to general N, though it requires ao- 
curate numerical integrals. Here we show that these methods could be useful tools for computing 
quantities in (11.11) and establish efficient numerical procedures. Our methods do not depend on 
specific distributions on N and Xi and therefore could be applicable under general settings. 

In the remainder of this section, we present a motivating application and its literature in Subsec- 
tion ll.ll and introduce notations used in Subsection ll.2[ In Section [2] random sums of discrete r.v.’s 
are treated, where computations of | 5 a?], k e N are investigated in Subsection 12. H and those 
of E\X\ I 5 A?] are studied in Subsection 12.2[ Both the recursion method and the Fourier method 
are investigated. In Section [3] we consider random sums of non-negative continuous r.v.’s, where 
we take the Fourier approach for computations of both E[N^ | 5 a? < Jc] and E\X\ | 5 a? < 2c], a > 0. 
Finally in Section 01 numerical examples are given, which show that proposed methods work rea¬ 
sonably. As applications, we consider predictors for both the Poisson shot noise process and the 
compound mixed Poisson process. 

1.1. Motivating application and its literature. A motivating example is prediction in the Pois¬ 
son shot noise process of the form 



( 1 . 2 ) 


t > 0, 


1=1 


where Q < < T 2 < • are points of a homogeneous Poisson N{t) with intensity d > 0 and (L,) is 

a sequence of iid Levy processes independent of (F,) and such that L,(t) = 0 a.s. t <0. The process 
of this type has many applications in rather different areas (see O, ||29]| [l2n and ifTSll ). One of 
important research topics is the prediction of future increments M(t, t-l-5] := M{t+s)-M{f), s,t > 0 
based on the present observation M(t). For example, in non-life insurance M(t, t -l- 5 ] is interpreted 
as the number or amount of future payment in the interval (t, t -l- 5 ] from an insurance company to 
the insured. Another interpretation is that M(t, t -l- 5] may describe the workload to be managed by 
a large computer network for sources in the interval (t, t -l- 5 ]. Due to the properties of both Levy 
and Poisson processes, the prediction of future increments M(t, t -l- 5 ] given M(t) reduces to 

(1.3) E[M(t, t + 5] I Mit)] = E[N(s)]E[Liit + s - U)] + E[LM]E[Nit) \ 

where f/ is a uniform r.v. on (0, t) denoted by f/(0, t) independent of (L,). The proof of (11.31) is 
given in Appendix lAI or IfT^ (2.1)]. Here computations of £’[A(5)], £’[Li(5)] and E[Li(t + s - U)] 
are trivial. Since points (T,) of Poisson have the order statistic property, we can regard the sequence 
(Ti) in the quantity E[N{t) \ M(t)] as that of iid U (0, t) r.v.’s. Accordingly, taking X, : = L,(t- T,) and 
N := N(t) of M(t), we obtain the form (11.11) . Similarly, higher conditional moments E[M(t, t + s]^ \ 
M(t)], k e N are obtained as functions of £'[A*^(0 | 

A series of papers Il20ll . IfTOl and IfT^ assumes particular marginal distributions for (L,) such 
as Poisson or negative binomial, and exploits their specific properties to obtain the conditional 
moments. Although asymptotic behaviors of E[N \ Sn = k], k ^ 00 have been studied in lfT0]| and 
||25][2^ . only limited distributions are treated. Our methods presented here require no particular 
assumptions on distributions of N and Xi and therefore could be applicable under more general 
settings than those of previous papers. 
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1 .2. Necessary notations and tools. Throughout we use the following notations related with gen¬ 
erating funetions. For a fixed r.v. X and a non-negative funetion / and \u\ < 1, 

CO 

Gxiu) := Gf(u) := G,f(u) := 

k=0 


I 


u^df(x). 


where the last one is defined as a Riemann-Stieltjes integral if exists. From these quantities we ean 
obtain the Fourier (-Stieltjes) transforms </>{.)(«) = G(.j(e"0- Note that we use generating funetions 
not only for r.v.’s but also for diserete sequences (see OTl f. though after a proper standardization, 
they are the same. We write No := {0,1,2,...} and R+ := [0, oo) in the sequel. Moreover, braces { I} 
denote the Stirling numbers of the second kind (see [[B p.824]): the number of ways of partitioning 
a set of n elements into k non-empty sbsets. 

We say that the probability mass function qn = P{N = n) belongs to the Panjer (a, b) class if it 
satisfies 




n ~ 


n' 


q„-i, n 6 


for a + b > Q and a < I ( 112411201 p.l22]). Poisson, negative binomial and binomial distributions 
belong to this class. For later use, we present the Panjer recursion formula (see [|24ll20ll for details 
and the proof). 


Theorem 1.1. Suppose that N belongs to the Panjer {a,b) class and denote an lid sequence of 
non-negative integer-valued r.v.’s by (X,). Then 

P{Sn = 0) = E[P(X^ = 0)"^], n = 0, 

= 1 - ap'x, = 0) B(" + " - *■ 

Here we let 0° = 1 conventionally. 


2. Random sums of discrete random variables 

2.1. Estimation of number of iid components. In this section calculations for conditional mo¬ 
ments E\N^ I S n], k e N will be investigated, where iid random components (X,) are integer¬ 
valued. In case N belongs to the Panjer class, we apply the recursion formula to the calculation of 
\ Sn]- For general N, we consider the generating function of 15^^ = -] and then apply 
the inversion formula. 

Throughout we denote the expectation of a r.v. X over a measurable subset A c Q by E\X-,A\ = 
E[Xfx^A)\ following Kallenberg [[111 p.49]. Since we obtain P{Sx) by Theorem I l.li we mainly 
consider k e N which yields EiN'^ \ = E{N^\SnMP{Sn)- 

Theorem 2.1. Let N belong to the Panjer {a,b) class and iid r.v.’s (X,) take values in Nq. Let 
Co := aP(Xi = 0). Then, the restricted moments mf,{C) := E{N^',Sxi = (^1, k, C G No satisfy the 
recursion, 

m,(0) = E[X'P(Xi = 0)""], 

mm = T^fCo £ (^) mm + Xi (« + = 2) Xi (^) ^ ^ 

where mo(T) = P{S n = £)■ 
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Since for the ealeulation of a eombination of i <k - \ and m,( 7 ), i < k, j < i - \ 
is suffieient, we can recursively ealeulate the quantity. 

Proof. The iidness of (X,) and independenee between N and (X,) yield 

mm = E{N^-,S^ = 0] = £[iV'£[/(5„=o) I A^]] = E{N'^P{X, = 0)^]. 

Next we consider mk{£), £ >\. Conditioning argument and the Panjer (a, b) elass assumption yield 


(2.2) mm = Yu ‘ = Z ‘ ^ 

i=\ i=i ' 

where = P{N = n), n e Nq. Sinee (X,) are iid r.v.’s 


i-l 


1 


{a + ^) = a + b-YE[Y I S t = £] =a + bE[^ \ S,- = £] = E[a + b^ I Si = i]. 




i=i 


-Si 


Moreover, 
bXi 


[a + I S; - ^ (a + -j)P{Xi - j I A; - 0 ^ (a + y) 

i=o j=o 


bj^P(, X,=j)P(,Si-i=£-j) 
P(Si = £) 


Substitution of this into {a + b/i) of (12.21) and multiple interchanges of the order of summations 
give 


oo i 




'■=1 j =0 
e 


Y(a + ^)P(Xi = 2) j E ^(^-1 = ^ - 2>Vi) 


y=o 


1=1 


=E ("+ =i)Y,F{s.=e- j) x; m .-v, 

y =0 ^ (=0 h=Q ^ ' 

= E E (^) 

i=n h=n V ' 


k-i 


/i=0 


= aP(Xi =0)m,((^) + aP(Xi =0)^(^)m,(O + E(« + y)^(^i 


y=i 


/i=0 


Thus we obtain the desired result. 


□ 


If we take k = I with ^ > 1 in Theorem 12.1[ a rather simple expression is obtained 
1 ^ b' 

mm = y^jCoP(S^ = 0) + (a + -j)P{X, = j){P{St, = £ - j) + mf£ - y))}, 

which together with Theorem 1 1.1[ yields the conditional expectation. 

Next we consider the generating function for mk{£) with N a general r.v. 
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Proposition 2.2. Let N be a r.v. on Nq and let (Xj) be an iid sequence of r.v.’s on Nq. Assume 
< oo, e N. Then the generating function of the truncated k-th moment m^^if) = S = 

has the form 

(2.3) G^,{u) = Z { y } 4.(“) G^^{Gxfu)\ \u\ < 1, 

where the quantities by braces {} denote the Striling number of the second kind, and Gy\u), 7 e N 
denotes the j-th derivative ofGyix) at x = u. G^j\u), 7 e N denotes the j-th derivative ofGyix) at 
X = u. 

Proof A direct calculation yields 

00 00 

GmM) = Z = ^] = Z = i\N) = EN’^Glfu), 

f=0 i=Q 

where EN^ < 00 assures Fubini’s theorem since \Gxfu)\ < 1. We use the relation of the falling 
factorial {x\ = x{x - 1) • • • (x - /: + 1) and 

(2.4) = { 0 } = ^’ 

namely, 

EVN’^Glfu)] = Z { y } § i ^ } G{^{u)E[{N)jGl-J{u)] = dH, 

where we change the order of derivatives and the summation, which is valid from EN^ < 00 and 
\Gxfu)\ <1. □ 

In order to obtain mkif) = E[N^-, Sn = from G^*, two methods are considered. One requires 
numerical integrations and the other needs derivatives of at the origin. Since \Gm^{e"’)^ < 
{EN^f < 00 , we have G,„fe'^) e L^{-n,n). Then the Fourier expansion of Gm^{e"’) is guaranteed 
and their coefficients satisfy formula 

(2.5) mm = ^ f _ e -‘^‘ Un ,( e ‘“) du , £ e No, 

which correspond to the inversion of the Fourier transform G^m^')- *^he other hand, if we 
take derivatives of G^^ at the origin, we obtain mk{£) = ^Gm^(O). In view of (12.31) . however, 
the calculation of G^] would yield additional complexities, though we may possibly find some 
efficient recursion methods for a limited class of N. The choice of the two methods depends on 
distributional assumptions on N and Xi and we need numerical experiments to judge which is 
better. 

2.2. Estimation of magnitude of each iid component. In this subsection we consider the ex¬ 
pected magnitude of r.v. X\, k 6 N under the observation of the total number S n- Since the condi¬ 
tional moments minimize mean squared errors, we will consider= £'[Xj | S^^], k e N. Since the 
direct application of the Panjer recursion seems difficult for;^^^ and easy forXk+ ■= E[X\ \ S a^+i], 
we derive the recursion only for;^f^+. Meanwhile, the Fourier approach is applied to both. 
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Theorem 2.3. Let N be a Panjer {a, b) class distribution and let (X,) be a sequence ofiid r.v.’s on 
Nq. Assume EX'^ < oo, k € N, then the truncated k-th moment Xk+iO = E[X\',Sn+i = ^], ^ 6 N 
has the form 

XkAD = P(Xi = l)E[PiX, = of], and for £ > 2, 

^-1 L -k 

Yj {<^Xk+f - j) + - 7))- 


Xkfi) = = i)P{SN = 0) + 


1 


1 - aP(Xi = 0) 


y=i 


- ] 


Proof. For i = I, due to the iidness of (X,), 

Xk^a) = E[X\-,Sm^i = 1] = F(Xi = l)E[P(X, = Of], 
Let Cl := 1/(1 - aP(Xi = 0)). For £ >2, the property of N yields 
Xkf£) = E[X\-,SN^i=£] 


= ^fP(Xi=j)PiSM = £-j) 


,/=i 


e-i 


t-j 


fPiX, = £)P{Sm = 0) + Yj Z = £-j-m) 


7=1 


m=\ 


i-i' 


e-\ t-j 


= ^'F(Xi = £)P(Sn = 0) + Cl] a ^ fP(X, = j)P{X, = m)P{S= £ - j - m) 


j=l m=\ 
e-\ .V 

J 


+ bj] Y^.PiXi = j) Z >^P(Xi = m)P{SN = £-j-m) 

7=1 ^ m=\ 

= fPiXi = £)P{Sn = 0) + Cl] a ^ PiXi = m)E[X\-, 5iv+i = ^ - m] 


m=l 


.k \ 

+ Z = 7)£[Xi; Siv+1 =£-j] [, 

j=\ J i 


where in the third step, we use the Panjer reeursion for P{Sn = £ - j)- Finally, we arrange two 
sums and obtain the result. □ 

For the ealeulation ofxk = E\X\ | S a?] a direet applieation of the Panjer reeursion seems diffieult 
and alternatively we try the Fourier methods. For this we need the generating funetion oixk- 

Proposition 2.4. Let N be a r.v. on No and let (X,) be an iid sequence of r.v.’s on Nq. Assume 
CXj < oo, k eN, then the generating function of the truncated k-th moment XkO) = C[Xj; Sat = •] 
has the form 

(2.6) G,,{u) = Z { ^ I 


,/=i 


GxXuy 


\u\ < 1. 


Xk{£) = E{X \-Siv = ^] = Z /^(^i = j)PiSN-i = ^ - j). 


7=1 


Proof. In view of 
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the functionis the convolution of two non-negative functions giU) '■= fP{^\ = j) and giij) '■ = 
= j). Since = i;y=! u^P{Sn-\ = j) = £[GJ^‘(m)] = G 5 „(m)/Gx,(m), the generating 

function of g\{j) is enough. We use the relation of the falling factorial (12.41) and obtain 

E[xy^] = I ^ = EI y 

where we apply Fubini’s theorem, which is possible by EX\ < oo. Now the product of G 5 ^_| and 
E[X^u^^] yields the result. □ 

Similarly as before, two methods are considered to obtain from G^^. One is to use deriva¬ 
tives at the origin,= Gfl{0)li\,i > 1. The other is the inversion of generating function 

Xk{i) = ^f_ € € No. 

In view of (12.61) . the former method requires some efficient algorithm for calculating derivatives of 
G^j, whereas for the second one, accurate numerical integrations are inevitable. 


3. Random sums of continuous random variables 

In this section, we assume continuous distributions for an iid random sequence (Xj) taking values 
on R+, while keeping N to be r.v. on Nq. Similarly as before we consider E[N^ | S and E\X\ \ 
5 A?], k e N. Here the Fourier Stieltjes transform (FST for short) is our main tool. 

3.1. Estimation of random number from random sum. We firstly consider E{N^ | Sa? s 
[0, a]] = E\N’' I Sa? < a] for a g R+ and k G N. We are starting to observe the integral equa¬ 
tion as in [l24l Sec. 4.4.3], which corresponds to the recursion formula when is a discrete 
distribution. 


Theorem 3.1. Let N be a Panjer {a, b) class distribution and assume iid r.v.’s (X,) take values on 
R+ with common distribution Fy,. Then the restricted k-th moment to the Borel set by {S ^ < x}, 
mk(x) = E[N''; Sat < a], k g N satisfies the integral equation, 

(3.1) = ainik * Exfiix) + j*Exfi(x), x>0, 

where the operation * denotes the convolution as usual. 


Proof. Since N belongs to the Panjer (a, b) class, we can write 

00 CO CO 

mfix) = 'Yjn’^qnPxfix) = a^{n + if qnP*j’'*^\x) + b^{n + lf~^qnP*ffJ^\x), 

n=0 n=0 n=0 

where P*f{x) denotes the distribution of the n-th convolution of Xi. Using the binomial expansion 
and changing the order of summations, we obtain 

oo k /I \ oo k—\ 


coA:/y\ oo K—i /J 

mk{x) = « E E ( •) ^^^nE*x*^\x) + ^ E E ( • I 

n =0 ;=0 n =0 ,/=0 \ J ' 

CO oo r k—\ /J \ I 

= a E + E T E ( /) ^ 


n=0 


n=0 ^ 7=1 


k-1 


Xi 

k-l 


,/=0 


J 


n-’q„P*f’^^^\x) 
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= a(mi * FxC){x) + ^ ^ j| ^ n^qnF 


( n + l ) 
^1 


{x). 


□ 


In view of expression (I3.1L the integral equation seems useless to obtain m^ix) and we need ad¬ 
ditional techniques such as discretization of the density function of Xi as in If24l Example (p.l23)]. 
However, it is helpful to obtain the generating function of by providing an efficient recursion. 

Lemma 3.2. Assume that N belongs to the Panjer (a, b) class and lid r.v.’s (Xi) take values on R+ 
with common ch.f Then the FST of m^lx) = Sat < x], k gN has the following form 

k-\ 


fmfF) — 


1 


1 


afxfu) 


E 


-I- b 


k-l 

j 


(t)mfu)(pxfu), U£ 


The proof of Lemma is a straightforward calculation and we omit it. Notice that due to Lemma 
13.2[ (pmjji) can be presented by a combination of Gx/ifxfu)) and (pxfu) since fmoM = = 

Gn(<PxM)). Lor a general N, we directly calculate the LST of m^.. 


Proposition 3.3. Let N be a r.v. on No and let (Xf) be an iid sequence of r.v.’s on 
EN^ < oo for k gM, then the FST of mk(x) := 5 a? < x] has the form 

k 


Assume 


(3.2) 


f 


e‘^^dmk(x) 




E 

,/=i 


0Zi(“) 0 v(0Xi(m)) 


where the left integral exists in the sense of the improper Riemann-Stieltjes integral. 

Proof. Observe that m^.(x) is a bounded non-decreasing function and e'^^ is continuous for every 
M G R, then a Riemann-Stieltjes integral e'“dm;.(x) exists for all M > 0 [1301 (2.24) Theorem]. 
Moreover, integration by parts (twice) and Lubini’s theorem yield 


J r*M 

„' 


’‘dmk(x) 




J r^M 

iue“ 

0 


’‘mk(x)dx 




[e'“m,(x)]o"-L 




J r^M 

iue'“ 

0 


'P(Sn < X I A^)djc 




F^’^PiSx < X I iV)]" 


= E 




J nM 

. ' 


J r^M 

e'“dP(5^ < X I A^)} 

0 


e'“dP(S^ < X I A^) 


where in the third step, we use m^ix) := E[N’"P{S at < x | A^^)] for all x > 0. To obtain the first 
equality of (13.2!) take the limit M —> oo on both side, where in the right-hand side, the limit and 
expectation are exchangeable due to EN^ < oo. Since the third equality follows similarly as in the 
proof of Proposition [2]2l we conclude the result. □ 

The inversion of LST is well known [fT2l Theorem 4.4.1]: Let f(u) be the LST of a bounded 
non-decreasing function F(x), then is defined as 

1 

F(x) = r~\(p(-) ](x) = lim — - 1)/ - /t]0(t)dt, x > 0. 

T—^oo '2,7T ^—Y 
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Hence, mjt(A) = T H ](t). When a r.v. is a Poisson with parameter A, so that a = 0, b = A 
in Theorem 13.11 we obtain mi (a) = A(mo * Fxi )(x) and m 2 ix) = A(mi * Fx^Xx) + mi(x). Thus, it 
follows that 

mXx) = ]W, m2ix) = ]W + m^x). 

3.2. Estimation of magnitude of each iid component. A direct application of the Panjer recur¬ 
sion seems difficult for both;^f^ and;^^i+, and we alternatively invert FST of these functions. In 
order to obtain the FST, we represent dindiXk+ in the form of a convolution. 


Lemma 3.4. Let N be a r.v. on No and (Xf) is an iid sequence of r.v.’s on R+ with common ch.f. (f>x^ 
such that EAj < oo. ThenXk(x) = E[X^^; Sx < a], e N has the form, 

xki.x) = 


Proof. We exploit the expression 


Xk{x) = E{X\-Sx<x\= f y^P(Sx-i<x-y)dPx,(y), 

Jo 

which is the convolution of P(Sn-\ < •) and £ y^dPxfy). Since the ch.f. of S is <psAF)l4>x{u), 
the conclusion is implied by the FST of £ y^dPxfx) which is e‘“^x^dPx,(x) = where 

EXf < oo assures the existence of fxl^u). □ 


The corresponding result for ElX’l; S< x] is obvious. Under the same condition of Lemma 
IMI we have;^-*+(x) = 


4. Numerical Examples 


We prepare notations of distributions used in examples. Denote a Poisson distribution with 
parameter A by Pois(A) and by Geo{p), a geometric distribution with parameter p of which proba¬ 
bility is P{X = k) = pq^, q = \ - p, k £ Nq. As usual write A ~ • if r.v. X follows the distribution 
after the tilde. All computations are done with Mathematica ver. 9 of Wolfram. 

Firstly a simple example of E[N | Siv] is presented by setting N ~ Pois(A) and Ai ~ Pois{y). 
We examine two proposed methods for mi = E\N\Sx], the recursion method and the Fourier 
inversion. For the probability of Sn, we use the ordinary recursion ('Theorem ll.il) . which yields 


(4.1) 


P{Sn = 0 


E[P(Xi = 0)^] = 

Z5=i "iP{x,=j)P{Sx = i-i), 


We apply Theorem 12. II to obtain the recursion. 


i = Q, 

e>\. 


(4.2) 


mfb) = 


E[NP{Xi = 0)^] = Ae-^eF<^^-^\ b = 0, 

Z5=i T^(^i = j){PiSN = i-j) + mfe - j)}, i>\. 


Another method for mfi) is to apply (12.51) to the Fourier transform (Proposition l2.2l) . which is 


G,„,(e'“) = ENGlfeA = 

In FigurelH we plot E{N | S a? = ("I = mfi)IP{S n = ^), ^ > 0 using both methods. Although they 
coincide when parameters are moderate, if either of parameters of Poisson for N and Ai is large, 
we observe instability for small b in the Fourier approach (Figure [U Right, squared dots), though 
for large £ there is no difference. 
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Figure 1. Left: the eonditional moment E[N \ S n = ^ & [0,400] when Poisson 

parameters for {N,X\) are (4 = 20, y = 10). Right: E[N | Sa? = ^], ^ e [0,1000] 
under the setting (d = 100, 7 = 5). Squared dots are values by the Fourier approach 
and round dots are those by the recursion method. In the former case values of both 
methods coincide. However, in the latter case instability is observed in the Fourier 
method for small though for large i they coincide. 


Next, we consider an example of E[Xi \ S /v+i] for the recursion and that of E[Xi \ Sn] with 
the Fourier transform, where N ~ Pois{A) and Xi ~ Geo{p). The Panjer recursion is applied 
to both P{S N = ■) and P(S n+i = ■), for the latter of which we also use the convolution. For 
Xi+(0 = 5AT+i = £], we use the recursion by Theorem [23] i.e. 

( pqe^'^‘^, ^ = 1, 

= I + Y!-\ - j), i > 2. 

Eorxm = E[X, , S — ^], the inversion of the Fourier transform (12.5! ) is applied to 


G,ny“) 


qe (^^(l-e"')i 

--exp d-^-— . 

1 - qe"^ ^ qe'^' - 1 > 


In Figure [21 we plot E\Xi | S a? = ^] and E\Xi j S a?+i = i] for £ > 1. Since the graphs show very 
similar curves for a moderate setting of parameters, we conclude that both methods work properly. 
However the instability is again observed in the Fourier approach (Figure [2l Right, squared dots) 
when the parameter d is large and i is small. 


4.1. Prediction in Poisson shot noise process. We pursue the prediction £’[M(t, t+ 5 ] \ t,s> 

0 of the model (II.2L i.e. calculate the quantity E{N{f) \ M{t)\ in (11.31) . As mentioned, since the 
order of {Tj) in (11.21) does not change the distributional relation of N{t) and M{f), by the order 
statistics property of the Poisson, we may consider the model M{t) := Gk{t - Uk) with the iid 
U(0,t) sequence (t/,), and then study E[N(t) \ We assume that the processes L^’s are iid 

compound Poisson processes such that the generic process L has the form L(t) = Yj^ where 
No(t) ~ Poisiyt), and (Ty) denotes an iid sequence of non-negative jump sizes. 

Now by setting N := N(t) and Xi := L,(t - f/,), i e N, the calculation of E[N(t) \ M(t)] can be 
considered in the framework of E[N j Sa^]. For the probability of Xi, since Ao(t - f/i) does not 
belong to the Panjer class, we take the Fourier approach. For this we need ch.f. of Ai ■= Li{t-U\), 
which is 


iuXi^ _ _ 


= E{e‘ 


^] = 


(«)-!) _ 7r(c^yj(H)-l) 


7(1 -0Fi(w)) 
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E(X, I = 0 A E(X, I = 0 




( 


Figure 2. Left: the eonditional moments E\Xi | S a? = ^] (square dots) and E\Xi \ 
S AT+i = £] (round dots) for £ e [1,300] when parameters of Pois(A) and Geo{p) for 
{N,X\) are (d = 40, p = 0.25) respeetively. Right: the same quantities of the left but 
with ^ e [1,1000] and (d = 150, p = 0.2). In both graphs these quantities present 
quite similar eurves for large £. However, in the right graph instability oceurs in 
small £ of E\Xi | S at = ^] (the Fourier approach). 




0 5000 10000 15000 0 50 100 150 200 250 

( ( 

Figures. Left: we plots £'[A^(1) | M(l) = £],££ [0,18000] which is the quantity in 
the predictor of a Poisson shot noise process of (11.31) with 1=1. Right: we plot the 
predictor of a compound mixed Poisson process, £’[Z(1,2] | Z(l) = £],£& [0,250]. 
In both graphs one see the non-linear curves which show that the linear estimations 
are insufficient. 


where (pYiiu) is the ch.f. of Y\. Thus after putting Ti ~ Poisip) so that (pYiiu) = we 

obtain the probability of Xi by the Fourier inversion. For simplicity, we set t = 1, i.e. consider 
E[N(l) I A/(l)], and apply Theorem |2.II or equivalently apply the recursions (14.11) and (14.21) with 
initial values P{Sn = 0) = e'HP(Ai=o)-i) mi(0) = E[NPiXi = 0)^] = dP(Zi = 0)P{Sn = 0). 

In Figure [3] (left), we plots E[N(1) \ M(l) = £] for £ G [0,18000] with y = 100, p = 5 and 
d = 30. In view of the graph, our computational method seems to work well, and one can see a 
non-linear curve which shows that the linear estimation of N{t) by M(t) is insufficient. 

4.2. Prediction in compound mixed Poisson process. We consider an example of the compound 
mixed Poisson process mixed by a Gamma r.v. called compound Polya process 0 Ex. 4.1]. Let 
N(t) := n(6A(t)) denotes a mixed Poisson process where n(t) be a homogeneous Poisson process 
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with intensity 1 on [0, oo), A(t) is an intensity measure and 0 is a Gamma (a,/3) r.v. of whieh 
density is fe^x) = . Then the proeess has the form Z{t) = t > Q, where 

Ay’s are iid r.v.’s on Nq or R+ sueh that N and {Xj) are independent. Since the cr-fields Qt by 
{N{t), N(t, t + 5 ], Z{t)}, t,s > 0 and by {N(t), Z(t)} are finer than that by {Z{t)}, the conditional 
expectation of increments Z(t, t + 5] := Z(t + s) - Z{t) given Z{t) has 

E[Zit, t + s]\ Z{t)] = E[E[E[Zit, t + s]\g,]\ m I Z(0] 

(4.3) = E[Xi]E[E[N(,t, t+s]\H]\ Z{t)] 

= E[X,]E[E[N(t,t+s]\m)]\Z(t)], 


where in the third step we use the conditional independence of N(t, t + 5 ] and Z(t) given N(t) (' IfTTl 
Prop. 6.6]). Since 


E[N{t, t + 5] I N{t) = m] = J^k 


k=0 


t + 5 ] £[6/^+»^e-eA(f+^)] 

~k\ £[5/»ig-0A(f)] 


= A(t, t + 5 ] 


= A(t, t + 5 ]- 


a + m 


where in the second step we exchange the infinite sum and the expectation operator (see also lEl 
(1.4)]), we proceed the calculation (14.31) to get 


E[Z{t, t + .] I Z(0] = + E[N{t) I Z(0]). 

A(t) +yS 

Now let A(t) := t, jS := 1, or := 7 and X^ ~ Geo (1 /4), we obtain 

E{z{t, t + 5 ] I z(t)] = ^(7 + Emt) I zm. 

Since N{t) does not belong to the Panjer class, we apply the Fourier approach. Due to Proposition 
l2.2l together with 


GsAu) = E{Glf{,u)] = 


1 

\+t{\-GxM)) 


and G,„,{u) = E{NGl{u)] 


7Gzi(m) 

{l+t(l-Gx,(i^))}«’ 


we obtain the quantity E\N{f) \ Z(t)] by the inversion formula (12.51) . In Figure [3] (right), we plot 
E[Z{t, t + 5 ] I Z{t) = £], i £ \0, 250] with s = t = 1, where one would again observe a non-linear 
curve. 


Appendix A. Calculation of (11.31) 

For the calculation of (II.3L we use the following properties. 

(a) . By definition, the cr-field by M{t) is included in the cr-filed by - T^)) and {Ti£). 

(b) . Since A(0 = Zr=iW)> the cr-field by N{t) is included in the cr-field by (T^). 

(c) . By the order statistics property of a Poisson, given N{t) and N{t + s), the set of points 
(r^.) G (0, t] and the set of points {T^) &{t,t + 5 ] are independent. 

(d) . Given N{t, t + 5 ], points (T^.) £{t,t + 5 ] are mutually independent. 

(e) . Stationary and independent increments of Levy processes. 
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By a multiple use of iterated property of the eonditional expeetation [fTTl Theorem 6.1 (vii)], de¬ 
tailed the ealeulation of (11.31) is 

t + 5] I M(t)] 

Nit+s) N(t) 

= E[Yj -Tj,t+s- Tj] I M(0] + Lj(t -Tj,t + s- Tj] \ M(t)] 

j=N(t)+l 7=1 

N{t+s) 

= E[ Yj E[Lj(t -Tj,t + s- Tj] I N(t), N(t + s), m), (Lu{t - n))}t.T,<t] I M(t)] 

j=Nit)+\ 

Nit) 

+ E[Y E[Lj(t - Tj, t + s-Tj]\ mi {Luit - n))m<t] I M(t)] 

J=i 

Nit+S) 

= E[ Y E[Lj{t + s - Tj)mj<t^s) I N{tl N{t + s)] | M{t)] 

j=Nit)+\ 

Nit) 

+ E[Y E[Lj(t -Tj,t + s- Tj] I Tj, Ljit - Tj)] \ M(t)] 

i=i 

Nit+s) Nit) 

= E[Y E[L{t + s-U)]\M(t)] + E[YE[L{s)]\M{t)] 

j=Nit)+l 7=1 

= E[Nit, t + 5] I M(t)]E[L(t + s-U)]+ E[Lis)]E[N(t) \ M(0], 

where in the second step, the properties (a) and (b) are used, and in the third step, we exploit (c) and 
(e) so that the conditional independence of (Lj(t -Tj,t + s- T j]) j-t<Tj<t+s and {(T,), (T,(t - r,))},: 7 ’,<f. 
In the fourth step we use (d) and (e). Finally since the quantity N(t, t -l- 5 ] is independent of the 
(T-field constructed by all available set before t, the conclusion holds. 

Acknowledgment: The author would like to thank Prof. Tomasz Rolski for fruitful discussion of 
the topic and hospitality when he visited the mathematical institute in the University of Wroclaw. 
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